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Abstract 

Previous results on fermion chirality-flipping four-point functions are ex- 
tended to non-abelian gauge theories. The problem is purely non-perturbative, 
and the analytical formalism used is based on the Schwinger-Dyson hierarchy. 
This is truncated, and the resulting equations are solved numerically by relax- 
ation techniques. Taking the large-iV limit in SU(N) theories simplifies the 
problem substantially, allowing the formulation of conjectures on the behavior 
of n-point and chirality-flipping fermion Green functions for general n. 

In this work we study the dynamical generation of fermion chirality nipping 
four-point functions within the framework of non-abelian gauge theories. These 
are generated in a purely non-perturbative manner, so in our analysis we have to 
use an adequate formalism like the one based on the Schwinger-Dyson (SD) hierar- 
chy. The resulting equations are analytically intractable and very challenging even 
numerically. However, a certain limit in the theory allows a drastic simplification. 

Our interest is focused on two objectives. One is to estimate the critical value 
of the gauge coupling, below which the four-point functions are equal to zero. We 
would like to study this in a case where four-fermion condensates develop on scales 
higher than the scale of two-fermion-condensate formation. This would make it 
consistent to treat the four-fermion condensate problem independently of the mass 
generation problem. Our other goal is to study the momentum dependence of the 
fermion four-point function. If such an object could generate a fermion or a gauge- 
boson mass, then two or 4 lines of the four-point function would be closed off into a 
loop or into another four-point function. It would therefore be interesting to study 
the details of the momentum dependence, such as the relative size of the four-point 
function when different pairs of momenta are large. 

In the present study we use a one gauge-boson exchange approximation, where 
one gauge boson can attach to any pair of the four legs. The present work constitutes 
a clear progress with respect to our previous study [2], since, apart from considering 
a general non-abelian group, it includes a treatment of non-linearities, it does not 
neglect terms proportional to external momenta while at the same time exploring the 
full available momentum space, and in a certain limit allows us to draw conclusions 
on the behavior of general n-point functions. 

1 Talk given at the 17th annual MRST meeting in Rochester (May 1995), based on work in 
collaboration with B. Holdom [1]. 



We focus our attention on fermion operators that are purely chirality changing of 
the form i>L^tpR' l pL^'i'R + h.c, where ip L = ^ 1 ~ 75 - > -0 and ipR = ^ 1+ ^ ip. We constrain 
our study to just one fermion flavor. The four independent operators which have 
this property and respect parity are 2 

$a^i)i>a^i) = i> L a^i) R i> L a^i> R + }i.c. 

^(^A>^A> + ^A a 7 5 ^A a 7 V) = i'L^Ri'L^^R + h.c. 

^AV'^AV^ = ^AV^W^AV^fl + h.c., (*) 

where A a are the generators of the non-abelian gauge group we wish to study. From 
now on we will denote all operators of the form ijjFifjipF'ip, where F, F' are matrices 
with possibly non-trivial spinor and color structure, by F <g> F 1 . 

The vacuum expectation values of these four operators receive exclusively non- 
perturbative contributions, and their non-zero values would break dynamically a 
chiral U(1)a, were the latter not broken by instantons. 

The SD formalism relevant to these four-fermion operators was discussed in 
Ref. [2] and led to an equation shown diagrammatically in Fig. 1. It results from a 
truncation of the SD hierarchy achieved by approximating the four-fermion+photon 
vertex appearing in the equation by a four-fermion vertex and a photon attached 
via a bare vertex on one of its legs. Such an approximation neglects, among others, 
non-linearities that would become important in the infra-red. We will return to 
them later. The equation has been symmetrized to include diagrams with gluons 
connecting all possible pairs of fermions, so the right-hand side is multiplied by a 
factor of 1/2. 

We are considering a four-point function associated with the Green function 
{0\T{ijj a ijji3ijjpijj T }\0) . In momentum space, the four-point function receiving exclu- 
sively non-perturbative contributions is 

+O c ° l + or P {\ a I aP ® \ a I pT + A a 7 ^ ® A a 7p T ) + O c T olor X a a^ ® \ a a^ (2) 

where the four scalar functions Os+p, Ot, O c g l °p and O^ lor depend on 6 variables, 
which are all the independent and Lorenz-invariant combinations of the external 
4-momenta pi,...,p4, corresponding to the fermions with spinor indices a,/3,p,r 
respectively. We wish to develop the SD equations for these scalar functions. 

The functional integral operators appearing below are defined exactly as in 
our previous work [2]. We work in the Landau gauge which is popular in studies of 
SD equations; the gauge boson propagator reads = p- — ^pr~) • 



2 We choose to work in Euclidean space, so there is no difference between upper and lower 
Lorentz indices. 



By combining the results of the abelian case [2] with considerations pertaining 
on the color structure in a SU (N) theory, we have 
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where T l [0] = T l {O l ) and (9* are the form factors with loop-dependent arguments 
corresponding to diagrams i = A, F [2]. 

We are faced with a system of four coupled 4-dimensional integral equations in- 
volving functions of 6 variables. As it stands, the problem is analytically intractable, 
and even a numerical solution proves to be beyond our present means. 

By taking the large- N limit, Eq.3 reduces to 
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It is clear that in the large- N limit, the two form factors Os+p,t decouple from the 
other two O s ° + °p T . As regards the first two, it is obvious that only Os+p is going 



be non-zero, since Or receives a negative contribution from the dominant diagrams 
A and B. It is also apparent that O s ° + ° P can be non-zero only if O T ° l0r is non-zero. 
Having a non-zero 0^ lor seems highly unlikely, however, since it receives a negative 
contribution from the dominant diagrams C and D involving itself. Even if it were 
able to be non-zero, it would require a critical coupling much larger than the one 
required for Os+p and it should therefore not present a problem. 

The form that this system of coupled integral equations takes in the large- N 
limit allows us to argue for the omission of contributions coming from operators 
proportional to external momenta. With regards to the operators with structure 
Pi<j^ v ®p^a vX and p^p^o^ ®PkPl& XT , which have so far been neglected, we can see 
that the contributions they receive from the dominant diagrams having the same 
functions in their vertices are negative, so we do not expect the corresponding four- 
point functions to be non-zero. Similar arguments make us neglect operators with 
structure 0\ = p^p^o^ ® 1 and C 2 = 1 ®p i iP v ^o ilv (and trivial color structure). 

We now write down the integral equation for Os+p according to the previous 
discussion. We have 

Nn c 1 / k o A k O B 

/n lya / j4 7 a L v ^S+P^l , ^S+P^i \ fr\ 

where q = p!+p 2 , = 8"" - (fc ~^ ) ^ ) ~ p ' r and the left-superscript of k Of + B P is a 
reminder that these functions depend on the loop momentum. 

We go in a reference frame where p[ = —p~2, and the form of the kernel of the 
integral equation indicates that Os+p is a function of five variables, i.e. Os+p = 
O s +p(p° 1 ,\p-'i\,pl\p-k\,q ). We also have k Oj +P = k Oj +P (k°, \k\,p 4 , \p 4 \,q°) and 
k O B s+P = fc Ol +P (rf,|lTi|,fc ,|A;U ). 

The structure of the equation indicates that a special solution to the above 
integral equation is factorizable, i.e. Os+p ~ B(p1, \p~i\, Q°)B(p1, \pk\, 0°)- Physically 
this is expected, since the momenta flowing in diagram A are independent from the 
momenta flowing in diagram B. Setting p\= p\= p° and = \p~l\ = \p\, we have 

B( P \ M; f) = ^ / W, |*|; 9 ») 7 ,^ r ^_ (6) 

where we have separated q° from the other two variables with a semi-column, since 
it does not mix with loop momenta and can thus be treated as a parameter. 

Note that for q° = this reduces to the linearized SD equation for the fermion 
two-point function. For non-zero q°, the problem reduces to the one of a three-point 
function. 

We have so far neglected non-linearities coming in our problem. These come 
from diagrams containing an odd number of four-point functions and a gluon in a 
way that makes them important only in the IR regime. They therefore provide our 
equation with an effective IR cut-off. Their non-linear nature gives our solution a 
natural scale. These diagrams are not easily resummable, so we multiply the kernel 
of our equation by the function A(q°) = fc2+A2 J 2 — , where A is the effective IR 



cut-off in our theory provided by the non-linearities, in an effort to mimic their 
effects, and we normalize B(0, 0; 0) = 1. Our full solution will thus be given by 

Os+p = (j^j B{p\, \pi\,q°)B(p4, \p~4\,q°), where r is a constant larger than unity 
because of the numerical factors appearing in the loops by which non-linearities are 
introduced in the problem. 

Moreover, we change our variables from p°, k°, \p\, and \k\, to p = ^p° 2 + \p\ 2 , 

k = ^k° 2 + \k\ 2 , <p = arctan (\p\/p°) and 4> = arctan (\k\/k°), so B(p° ,\p\; q°) = 
B(p, (ft; q°). We then solve for the function pB in order to increase the accuracy of 
the numerical solution that follows, since we expect B to decrease with increasing p. 
After some Dirac algebra, and omitting terms with structure PiP^cr^ (g> 1 according 
to our previous discussion, the equation becomes 
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The form of the integral equation allows us to use the same discretization lattice 
for the arguments of the function B inside and outside the integral. We choose the 
number of lattice points in each dimension to be equal to 40. Our results do not 
change by varying this number. The integral equation is solved via a relaxation 
method. We first insert an initial configuration for our function, and then iterate 
the equation until it is satisfied to a reasonable accuracy. 

We start by setting q° = 0. This gives us a critical coupling Na c = 2.7±0.3. This 
coupling is somewhat larger than the one required for a non-zero two-point function, 
which, for large-A^, is given by the relation Na c = ^ « 2.1. This discrepancy is 
expected to disappear in the limit Auv — ► oo. We then keep this coupling fixed, 
and for different values of q° we compute the corresponding values of -6(0, 0; q°). In 
Fig. 2 we plot the function -| B(p, (j); 0). 

We see that j^B(p, <fi; 0) exhibits the cos (7 log (p/A)) behavior ( with 7 a coupling- 
dependent constant and A the IR cut-off of the equation) that is well known in 
two-point function studies, and that it is independent from </>, as expected. We do 
not trust the behavior of the function for momenta smaller than A, since we expect 
our formalism to break down at that region. As q° increases and gets larger than 
AB(0, 0; q°), it plays the role the non-linearities played before and acts like an effec- 
tive IR cut-off, so the shape of B is the same as the one of the two-point function 
having q° as its IR cut-off. The dependence on <p remains very weak. This behavior 
was expected for q° = 0, and it persists for larger q°. Moreover, B falls rapidly with 
increasing q°, so the non-linearities become increasingly irrelevant in this regime. 
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Figure 1: The schematic form of the SD equation. We have labeled the four fermions 
by their spinor indices. We label the diagrams by the capital letters A,...,F, and 
we omit the factor of 1/2 multiplying the right-hand side. 




Figure 2: The function ^ B(p, </>; 0) 



The behavior found is approximately B(0, 0; q°) & (l — jj^x) 2 • It is obvious that, 
for a given critical coupling, q° cannot exceed A by much before the whole solution 
collapses to zero. 

The form that the integral equation takes for large N suggests a possible gen- 
eralization of the above results. In particular, the diagrams that dominate are the 
ones in which a gluon is attached between two fermion fields with spinor indices 
contracted with each other, since color flows through a closed loop. Therefore, one 
might expect only condensates of the form | < (ipip . . .tpip + ifj r y 5 'ifj...'ifj'-y 5 i(j) >=< 
TpLipR-'-TpLiiR + h.c. >, where the ellipsis stands for n — 4 fermion fields paired with 
each other as in the first and last pair, to get non-zero expectation values at the 
same critical value of the gauge coupling as the 2- and four-point functions. 

Moreover, we expect the corresponding scalar n-point function Og +P to factorize 
and to be given by Og +P = A^.B 1 ...B k , where A c is the characteristic energy of the 
theory, and p = 4 — 3k with k = n/2. In analogy with the previous discussion, 
we have defined the functions Bi = B(pi, (fif, qi),i =,l...,k, where, as before, pi 
is the four-momentum of one of the fermions in each pair and qi is the total four- 
momentum of the pair. In order for each of the functions Bi to satisfy the same final 
integral equation as before, however, we have to assume that the three-momenta of 
the fermions in each pair are equal and opposite to each other, so that q} = 0. 
The qualitative results should not change substantially, nevertheless, if one makes 
the weaker assumption that we can go in an optimum reference frame where no 
fermion has a three-momentum which is much larger than the three-momentum of 
the fermion it is contracted with. 

To conclude, in this work we have treated a truncated Schwinger-Dyson equation 
for the four-point function of fermions interacting via a non-abelian gauge interaction 
in the Landau gauge. We considered chirality-changing four-point functions which 
receive non-perturbative contributions exclusively. The large- N limit of SU(N) 
renders our problem analogous to a three-point function one. Our numerical results 
indicate that the equations exhibit a critical behavior, and that the critical coupling 
appears to be roughly equal to the one required for the formation of two-fermion 
condensates. The beauty of the large- N limit is that it also allows us to guess the 
critical and momentum behavior of fermion chirality-flipping n-point functions for 
higher n. 
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